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The hadron energy resolution of a highly granular CALICE analogue scintillator-steel 
hadronic calorimeter was studied using pion test beam data. The stochastic term con- 
tribution to the energy resolution was estimated to be ~58%/^£;/GeV. To improve 
an energy resolution, local and global software compensation techniques were devel- 
oped which exploit an unprecedented granularity of the calorimeter and are based on 
event-by-event analysis of the energy density spectra. The application of either local or 
global software compensation technique results in reducing of stochastic term contri- 
bution down to ~45%/ ^ E/GeV . The achieved improvement of single particle energy 
resolution for pious is ~20% in the energy range from 10 to 80 GeV. 

1 Introduction 

To evaluate detector technologies for the future linear collider experiments [1] the CAL- 
^ \ ICE collaboration [2] has constructed and commissioned high-granular electromagnetic and 

hadronic calorimeter prototypes that have been successfully operated in test beam experi- 
^ I ments at DESY, CERN and Fcrmilab since 2006 till 2011. During the test beam campaigns 

a large amount of data was accumulated using muon, electron and pion beams in the energy 
range from 1 to 180 GeV with different detector configurations. The unprecedented gran- 
ularity of the CALICE calorimeter prototypes allows to analyze hadronic shower structure 
with high spatial resolution [3] and test particle flow algorithms using test beam data [1] . 

The hadronic energy resolution of the CALICE analogue scintillator-steel hadronic calo- 
rimeter was studied using test beam data. We also present a brief description of developed 
software compensation techniques as well as the results of their successful application to 
improve a single particle energy resolution for pions. 

2 Experimental setup and data 
2.1 CALICE test beam setup 

The data collected during 2007 test beam campaign at CERN SPS with positive and negative 
pion beams in the energy range from 10 to 80 GeV were analyzed. During this data taking 
period, CALICE test beam setup consisted of silicon-tungsten electromagnetic calorimeter 
(ECAL) 0, scintillator-steel analogue hadronic calorimeter (AHCAL) [5] and scintillator- 
steel tail catcher and muon tracker (TCMT) [7]. The test beam setup was also equipped 
with various trigger and beam monitoring devices including Ccrcnkov counter. 

The Si-W ECAL is a sampling calorimeter with a total depth of 24 radiation lengths. In 
this study the ECAL was used for event selection and early shower detection. Events with 
a primary inelastic interaction in the ECAL were rejected. 

The AHCAL consists of small 5-mm thick plastic scintillator tiles with individual read- 
out by silicon photomultiplicrs (SiPMs). The tiles are assembled in 38 layers (each layer 
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90x90 cm^) interleaved by 2-cm thick stainless steel absorber plates. The size of the scin- 
tillator tiles ranges from 3x3 cm^ in the central region, 6x6 cm^ in the outer region and 
12x12 cm^ along the perimeter of each layer. In the last eight layers only 6x6 cm^ and 
12x12 cm^ tiles are used. In total, the CALICE AHCAL has 7608 scintillator cells and 
amounts to a depth of 4.5 nuclear interaction lengths. 

The TCMT consists of 16 readout layers assembled from scintillator strips read out by 
SiPMs between steel absorber plates and has two sections with different sampling fractions. 
In this study the information from TCMT was used for muon separation and to minimize 
the effect of leakage from the AHCAL at higher energies. 

2.2 Calibration and event selection 

The calorimeter calibration procedures are described in detail in [5l|8]. The visible signal 
in cells produced by a minimum-ionizing particle (MIP) was studied with muons and is 
taken as a base unit of visible energy measurement. To reject noise, only signals above the 
threshold of 0.5 MIP arc used in the analysis, hereinafter these signals are called hits. 

The total deposited energy is obtained at the "electromagnetic scale" . The conversion 
factors from the visible signal in units of MIP to the total energy in units of GeV are 
extracted from electromagnetic calibration of each subdetector. The conversion factor from 
visible signal to deposited energy in the ECAL for non-showering hadrons was estimated 
using simulated muons and the measured muon response from test beam data sets. For the 
conversion to the hadronic energy scale, the ^ factor was considered in addition to calculate 
the initial reconstructed energy, as discussed below. 

The event selection procedure provides a purification of the pion samples that have an 
admixture of muons as well as electrons or protons. The efficiency of muon identification 
was checked using both simulated muon samples and muon samples from test beam data 
and was estimated to be better than 98%. The muon identification results in a purity of pion 
samples from muons better than 0.5% in the studied energy range. The Cerenkov counter 
is used to remove electrons from the tt" samples and protons from the 7r+ samples. 

To analyze the intrinsic AHCAL energy resolution for hadrons, an additional constraint 
is applied to the purified pion samples: a shower start (a position of the primary inelastic 
interaction) was required to be in the first five layers of the AHCAL. This requirement allows 
the effect of leakage into the TCMT to be reduced by selecting hadronic showers which are 
mostly contained in the AHCAL. 

2.3 Hadron energy reconstruction and intrinsic energy resolution 

Without software compensation the initial reconstructed energy -Einitiai of an event is calcu- 
lated as follows: 

-^'initial = -^'eCAL + ~ ' I -^hit + -E-TCMT | 7 (1) 

\m J 

where E^c^^i^ is the energy lost by a hadron track in the ECAL, |- = 1.19 is a scaling factor 
to take into account the different response to electrons and hadrons in the non-compensating 
AHCAL (the coefficient was obtained for pious by averaging the ratio of beam energy to the 
total energy reconstructed at electromagnetic scale over the studied energy range), i?hit is 
the energy deposited in one AHCAL cell and Stcmt is the energy deposited in the TCMT. 
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Figure 1: Reconstructed energy distributions for pious at 10 GeV (a), 40 GeV (b) 
and 80 GeV (c) without compensation (black circles and solid line) and after local 
(LC - blue triangles and dashed line) and global (GC - red squares and dotted line) 
software compensation applied. The curves show a Gaussian fit to the corresponding 
distributions. Statistical errors are shown. 



The reconstructed energy distribu- 
tions were fitted with a Gaussian in the in- 
terval of ±2 RMS around the mean value. 
Hereinafter, the mean and sigma of this 
Gaussian fit at a given beam energy are 
referred as a mean reconstructed energy 
-Eicco and resolution CTicco, respectively. 
Fig. [T] shows three examples of the energy 
distributions for three different beam en- 
ergies. The initial reconstructed energy 
distribution is shown by the black circles 
and its Gaussian fit is indicated by the 
black solid curve. 

The response of the calorimeter setup 
to pious as a function of beam energy is 
shown in Fig. [2ji, where black circles cor- 
respond to the mean initial reconstructed 
energy without software compensation. 
Due to the intrinsic non-compensation of 
the CALICE AHCAL, the response to pi- 
ous is non-linear with energy, deviating 
±2% from a perfectly linear behavior in 
the studied energy range. In Fig. [2]3 the 
relative residuals to the true beam energy 
are shown. The green band in Fig. [5] cor- 
responds to the systematic uncertainty of 
the reconstructed energy with the main 
contribution from the uncertainty of MIP 
to GeV conversion coefficient for the AH- 
CAL which was extracted from electro- 
magnetic calibration with an accuracy of 0. 
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Figure 2: (a) Linearity of response to pious 
and (b) relative residuals to beam energy ver- 
sus beam energy without compensation (black 
circles) and after local (LC - blue triangles) 
and global (GC - red squares) compensation. 
Filled and open markers indicate tt~ and 7r+, 
respectively. Dotted lines correspond to -Ercco = 
i?bcam and green band shows systematic uncer- 
tainties for the initial tt~ data sample. 
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The relative energy resolution is shown 
in Fig. [3J where black circles correspond to 
the initial energy resolution without soft- 
ware compensation. The resolution for 7r~ 
events (filled markers) is in good agreement 
with that observed for tt"*" events (open 
markers). The solid curve represents a fit 
with the following function: 

where E is in GeV, and a, b and c are 
stochastic, constant and noise contribu- 
tions, respectively. The noise term is fixed 
at c = 0.18 GeV, corresponding to the mea- 
sured noise contribution in the full CALICE 
setup taking in account contributions from 
ECAL (0.004 GeV), AHCAL (0.06 GeV) 
and TCMT (0.17 GeV). These estimates 
were obtained using dedicated runs without 
beam particles as well as using random trig- 
ger events constantly recorded during data taking [8]. 

The stochastic term contribution to the initial hadron energy resolution of the AHCAL 
was estimated to be ~ , The constant term contribution is ~ 1.6%. The similar 

^/E/GcV 

relative energy resolution was observed for Tr~^ and 7r~ samples. The application of soft- 
ware compensation methods described below allows to decrease the contribution from the 
stochastic term and to improve the hadron energy resolution for the AHCAL. 

3 Software compensation 

A response of sampling calorimeter to hadrons is much more complicated than to electrons 
and comprises the contributions from two different components: electromagnetic (mostly 
from tt" particles and nuclear photons) and hadronic. The latter includes a so called "invis- 
ible" component (neutrons, nuclear binding energy losses, etc.) that remains undetectable 
leading to lower calorimeter response to hadrons than to electrons. The fractional contain- 
ment of the mentioned components fluctuates significantly from event to event [QlllOj. The 
existence of the "invisible" component and event-by-event fluctuations of electromagnetic 
and hadronic fractions inside hadronic showers are two factors that deteriorate a hadronic 
energy resolution compared to electromagnetic one. 

For intrinsically non-compensating calorimeters, compensation can be achieved by using 
an approach known as "off-line" or "software" compensation and based on the expectation 
of higher energy density inside electromagnetic component compared to a hadronic one. 
Such procedures have been successfully implemented for several calorimeters, e.g. WAl 
calorimeter [11] , LAr calorimeter of HI experiment [12] and ATLAS hadronic endcap calo- 
rimeter [13] . 

To improve the energy resolution, two software compensation techniques based on hit 
spectrum analysis were developed for the AHCAL which allow to decrease the contribution 
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Figure 3: Relative resolution versus beam en- 
ergy without compensation (black circles and 
solid line) and after local (LC - blue triangles 
and dashed line) and global (GC - red squares 
and dotted line) compensation. The curves 
show flts using eq. [21 
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from stochastic term. Both approaches exploit a high granularity of the CALICE hadronic 
calorimeter and are based on event- by- event analysis of the energy density spectrum ob- 
tained from the signals of individual calorimeter cells. At the same time the approaches 
follow different ways to construct and apply compensation factors: in the local approach 
there are several weights applied to the signals of individual calorimeter cells, while in the 
global approach one compensation factor is applied to the energy sum calculated over all 
calorimeter. 

3.1 Local software compensation technique 

In the local software compensation procedure different weights are applied to every cell of 
the AHCAL. The values of these weights reflect whether a cell belongs to a hadronic or 
electromagnetic sub-shower. To distinguish between different types of energy deposit on 
the cell level, the local energy density in each particular cell is analyzed. As the cells in 
electromagnetic sub-showers typically have a higher energy deposit, they get a lower weight 
in the overall energy sum compare to cells with a lower energy density. The energy of the 
event corrected by the local software compensation technique can be calculated as follows: 

EhC = -E-ECAL + -^hit • '^hit(^^initial) H ' -EtcMT, (3) 

^ — ' TT 
hit 

where Whit(^'initiai) is a hit weight which depends on the hit energy density and on the total 
initial event energy. The procedure of adjusting weights and extracting parameters of their 
energy dependence is described in more detail in |14) . 

3.2 Global software compensation technique 

The global software compensation procedure uses one observable that describes the shape 
of the hit energy spectrum for a given event. For a higher electromagnetic fraction one 
can expect on average a higher energy deposition per calorimeter cell and therefore a larger 
relative contribution of high energy hits to the event hit spectrum. For each event the global 
compensation factor Cgi is calculated as a ratio of the number of shower hits with measured 
visible signal in MIP below threshold eum to the number of shower hits with measured visible 
signal below the mean hit energy in the given spectrum. The threshold enm is set to 5 MIP. 
The energy of the event corrected by the global software compensation technique can be 
calculated as follows: 

^^GC = -E'eCAL + -P(Cgh ^initial) ' I ^ -^hit + ^^TCMT 1 , (4) 

\ hit / 

where P(Cgi, -Einitiai) is a polynomial function to take into account the energy dependence of 
the compensation factor. Three parameters of this polynomial arc extracted from test beam 
data. The detailed description of the global compensation technique is presented in |15] . 

3.3 Application of software compensation to test beam data 

Both software compensation techniques described above do not require a prior knowledge 
of particle energy for the compensation to be applied. To take into account the energy 
dependence of compensation factors and weights, initial reconstructed energy is used for 
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the calculation of corrected energy on event-by-event basis. The parameters of the energy 
dependence of compensation factors or weights are derived using test beam data. To assure 
the sample independence of the methods, the selected test beam data set was split into two 
subsets: the parameters for both techniques were adjusted using one subset and then applied 
to another subset. It should be emphasized that the application of either local or global 
compensation approach does not distort a Gaussian form of the initial energy distributions. 

In fig. [5] the response to pious is shown after applying local or global compensation tech- 
nique comparing to initial response without compensation. In general, the response remains 
linear after compensation within ±1.5% that is slightly better than before compensation. 

The relative energy resolution before and after compensation is shown in fig. [3l A good 
agreement between the tt~ and 7r+ samples was observed. The application of software 
compensation results in a decrease of the stochastic term contribution while the constant 
terms before and after compensation agree within errors. Both compensation techniques 
show very similar performance, with the local software compensation providing a slightly 
smaller stochastic term. 

Fig. |4] shows a relative improvement of 
the energy resolution achieved by both de- 
veloped techniques. We define the relative 
improvement as the ratio of the resolution 
after software compensation crgc (local or 
global) to the initial resolution CTinitiai- As 
follows from Fig.|4l the improvement due to 
software compensation for test beam data 
ranges from ~12% to ^--^25% in the stud- 
ied energy range, the local approach demon- 
strating 3% higher improvement from 30 to 
60 GeV. The worsening of the relative im- 
provement at higher energies can be caused 
by a distortion of hit energy spectra due to 
increasing leakage into the TCMT. In the lo- 
cal compensation procedure, the shower hits 
from the TCMT are not weighted. In the 
global compensation procedure, the partial 
absence of hits in the hit spectrum results 
in wrong (or shifted) estimate of the global 
compensation factor. 
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Figure 4: Energy dependence of the relative 
improvement of resolution for local (LC - blue 
triangles) and global (GC - red squares) soft- 
ware compensation techniques. 



4 Conclusion 

The hadronic energy resolution of the high-granular CALICE hadronic calorimeter was 
studied using test beam data collected in 2007 on CERN SPS. The calorimeter comprises 
7608 cells in 1 m^ and is highly segmented in both longitudinal and transverse direction. 
The intrinsic relative energy resolution of the CALICE AHCAL for pious is ^ . 

y/E/GeV 

The contribution of the constant term was estimated to be about 1.6%. The unprecedented 
granularity of the CALICE AHCAL allows to apply a software compensation approach based 
on hit spectrum analysis to improve a single particle energy resolution. Two techniques were 
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developed and tested on data samples. The local software compensation technique applies 
different weights to individual cells while the global technique uses one compensation factor 
applied to the sum of cell signals. The achieved relative improvement of hadronic energy 
resolution varies from 12% to 25% in the studied energy range resulting in the reduction 
of the stochastic term contribution down to ^ , Both local and global software 

y/E/GcV ^ 

compensation techniques give almost similar gain in resolution. 

The proposed techniques do not require a knowledge of particle energy for the software 
compensation to be applied, the initial reconstructed energy is used as a first approxima- 
tion to calculate compensation factors and weights, which are energy dependent. Although 
this energy dependence places some restrictions on the implementation of both proposed 
approaches in jet reconstruction procedure, their application in the frame of particle flow 
algorithms is possible when disentangled single particle clusters are considered. 
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